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Volume III reports on tasks which are to be continued. Therefore, 
only a progress report for the fourth quarter is provided. Exceptions 
to this are tasks 2.8 and 2.12 which were terminated some months ago. 

A brief summary of the work under these tasks is provided here for 
completeness. 
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2.3 Field Measurements for Remote Sensing of Wheat 

Field Measurements activities during this quarter have included: 
project leadership and coordination, data processing, and evaluation, 
data library management, data analysis, and preparation for summer 
data collection at Williston, North Dakota. 

Project Leadership and Coordination 

Early in the quarter three sets of examples of field measurements 

I 

data were prepared. Examples of each of the primary data types were 
Included, with emphasis on the data from the helicopter and truck 
spectrometer systems. Two of the notebooks were forwarded to NASA/ JSC 
for NASA Headquarters and NASA/ JSC. On May 5 a presentation describing 
the project with illustrations of typical 'data, particularly spectro- 
meter data, was made at the LACIE Project Monthly Review Meeting. 

A plan and procedures for the sensor comparison experiments to be 
performed in Finney County, Kansas during mid-May were prepared and 
sent to NASA/JSC for review and comments. Inputs have been received 
from NASA/JSC that make the procedures more compatible for the instru- 
ment systems involved. The systems involved are the NASA/ JSC FSS 
and FSAS spectrometer systems and the Purdue/Lars Exotech 20C spectro- 

I 

meter system. The purpose of the sensor comparison experiments is to 
determine how well the output of the systems correlate. 

I 

Data Processing 

During this quarter significant progress was made in completing 
the processing of the 1975-76 Field Measurements Project data. The 
1974-75 data processing was completed the 'previous quarter. Processing 
of the fall 1976 data is currently underway. The status of the data is 


summarized in Tables 2.3-1. 
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FSS Data . The 1975-76 FSS data (both field averaged and single 
scan formats) were calibrated using the calibration panel spectral 
reflectance measurements determined by the truck-mounted spectrometers 
for the summer. An evaluation of the 1976 FSS spectral reflectance 
calibration tables and the 1975 FSS calibration tables prepared last 
fall indicated that the calibration of the 1975 FSS data could be im- 
proved upon. The 1975 FSS data were therefore recalibrated using up- 
dated tables. The change in the calibration of the 1975 data does not 
affect the relative comparison of spectra collected on a given data. 

The change, however, will affect relative comparisons from date to date 
in some cases. 

The processing of the fall 1976 data began during this quarter. 

These initial steps include key punching the ground observations. 

LARS Spectroradiometer Data . Processing of the Exotech 20C data 
collected at the Willis ton Agriculture E^eriment Station during the 
summer of 1976 was completed. A review of the data indicated that a 
few runs were of questionable quality and’ need to be reprocessed. The 
data will be distributed when the reprocessing of the questionable 
runs is completed. Processing of the Exotech 20C data collected in 
the "modeling field" is one- third completed. 

During this quarter several ways to reduce the amount of time re- 
quired to process the Exotech 20C data were examined. Some of the possi- 
ble ways include an onboard digitizer in the truck, minicomputer pro- 
cessing, and better procedures for handling the bookkeeping and ground 
observation data. At present the cost of the digitizer and minicomputer 
systems are restrictive; more evaluation is needed to determine the 
cost/benefits of the systems. New procedures for handling the book- 
keeping data and ground observations are being implemented this summer. 
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JSC FSAS Data . All of the spring-summer 1976 FSAS data has been 
received and processed. A review of the post-calibration of the FSAS 
field standard indicated that the calibration table for the 1976 FSAS 
field standard could be Improved upon in the 2.1 to 2.4 urn range. 
Therefore, that data is currently being recalibrated. 

Data Library 

The 1976 FSS field average and single scan format data were distri- 
buted to researchers at the Environmental Research Institute of Michigan 
and Texas A&M University along with the recalibrated 1975 FSS data. A 
report discussing the FSS calibration was distributed with the data. 

Also distributed to ERIM and Texas A&M were the 1975-76 crop year set 
of ASCS periodic observations and ground ‘truth inventories for the three 
test sites. Updating of Volumes 2 and 3 of the Field Measurements Data 
Catalogs for the 1975-76 and 1976-77 crop years is in process to re- 
flect the current status of the data. 

Data Analysis 

Analysis of the truck-mounted spectrometer data from the Garden 
City, Kansas and Willis ton. North Dakota 'Experiment Stations is well 
underway with data acquired in 1975 and has begun for 1976 data. Analy- 
sis of variance and covariance analysis are being used as tools to 
investigate treatment differences, evaluate crop separability, and com- 
pare the discriminability capabilities of the Landsat and proposed 
thematic mapper bands. Currently, an analysis is being carried out 
to quantify the change in reflectance of measurements taken at differ- 
ent times during the day; data acquired over a sample of test plots 
at the Williston, North Dakota site on July 16, 1976 is being used for 
this evaluation which is vital for the optimum design of the 1977 
data acquisition. 
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Graphs of the S-191 data from both the Fiimey County, Kansas and 
Williams County, South Dakota intensive test sites for 1975 have been 
examined to assess crop discriminability and changes in the spectral 
response of a crop over the growing season; statistical analyses will 
be performed to quantify these results. Single scan data are being 
analyzed to evaluate within field variability and compare it to between 
field and between crop variability. The same analysis tasks will be 
continued in the next quarter as well as regression analyses to corre- 
late spectral and agronomic measurements . 

In addition, the analysis of data set acquired over the modeling 
field at Williston, North Dakota on July 20, 1975 has been completed 
this quarter. The data set consists of spectral measurements over the 
range of 0.4 to 2.5 micrometers made at approximately half-hour inter- 
vals at five zenith and eight azimuth angles. A separate report of the 
results of this investigation is currently being prepared. 

Preparation for Williston, N. Dakota Data Collection 

Experiments were outlined and plans for data collection were devel- 
oped in preparation for this summers data collection activities in 
Williston, N. Dakota. The experiments include a continuation of spring 
wheat and small grain treatment plots and the modeling field activities . 
The past two years experience and analysis of the data Indicated that a 
change in the plot layout would improve the experimental design. The 
improved design is being Implemented for this summer. 

Also, during this quarter the Exotech 20C instrument was thoroughly 
tested. A new circular- variable filterwheel was installed for the sili- 


con detector. 



Table 2.3-1. Data Processing/Reformatting Status. 



1974-75 Data 

1975- 

■76 Data 


1976- 

77 Data 


Instrument/Data Type 

Completed 
& In Library 

Completed 
& In Library 

1 

1 

In 

Processing 

At 

JSC 

Completed 
& In Library 

In 

Processing 

At 

JSC 

Landsat MSS 

Whole Frame CCT (Frames) 

20 Frames 

61 

0 

H.A. 

14 

0 

N.A. 

Aircraft Scanner 
(Dates /Runs) 

19/149 

16/97 

0 

1 

2/10 

1/ 

2/ 

Helicopter Mounted Field 
Spectrometer (Dates/Runs) 

Field .Averages 

19/2,251 

27/2,193. 

0 


0 

4/ 

1/ 

Individual Scans 

19/35,000 

27,40,000 

0 


0 

4/ 

1/ 

Truck Mounted Field 
Spectrometer (Dates /Runs) 

FSAS 

6/65 

23/322 

0 

1 

0 

0 

3/ 

Exotech 20C 

24/1,577 

14/1,307 

5/1,361 


- 

- 

- 

Exotech 20D 

45/645 

•• 
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2.4 Scanner System Parameter Study 

Background 

The scanner parameter study was initiated in CY76 to conduct research 
on analytical methods of multispectral sensor system design. It was planned 
that basic information theoretic approaches would be taken to determine op- 
timum performance levels achievable in a given environment against which to 
test performance of actual or modeled systems. Scanner system modeling tech- 
niques were to be developed which would permit explicit evaluation of scanner 
performance without the use of simulation techniques. The study includes 
representation of the scene and information extraction process to provide 
models of the environment in which the scanner would operate. Thus three 
main activity areas were defined (scene, scanner, and information extraction 
modeling) and research was begun in CY76 

Shortly into CY76, the direction of the study was changed on request 
of the sponsor. A conventional simulation study of certain proposed thema- 
tic mapper parameters was requested using NASA 24 channel scanner data as 
input. This activity required almost all the project resources through the 
end of CY76 and the results were reported in the June 1976 final report and 
in LARS Information Note 110976 included as an appendix in this report. Re- 
search on the analytical model program was resumed in CY77. Progress has 
been made in the three basic areas defined for the study. The work is re- 
ported with respect to the appropriate tasks in the following sections. 

Task 2; Spectral Data Base Development 


The activities included in this task relate to the acquisition, re- 
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formatting preprocessing, error checking cataloguing and averaging of 
field spectrometer data for use in defining scene models. A large enough 
number of field spectra had been processed by the end of the third quarter 
to adequately support the planned research. The task was thus terminated 
at the end of the third quarter and was reported on at that time. 

Task 3: Statistical Model Definition 


The statistical modeling tasks have evolved as part of Tasks 4 and 7 

X 

and are reported as part of those tasks. This modification in task defi- 
nition was discussed at the end of the third quarter and Task 3 was thus 
defined as completed at that time. 

Task 4; Optimum Basis Function Study 

The selection of wavelength band sampling functions to maximize the 
performance of a scanner system for classifying certain scene classes is 
subject of this task. The wavelength sampling functions are being called 

J 

the basis functions here since the approach being taken is to represent 
the spectral response from the scene by a finite weighted sum of functions 
which are derived from the statistics of the spectra themselves. Two ap- 
proaches are being pursued. The first is based on the use of the Karhunen- 
Loeve (K-L) representation and the other is an information theory approach 
which utilizes the concept of mutual information. Both approaches have been 
discussed in detail in previous reports. A summary of progress on the two 
approaches in the fourth quarter is presented below: 

K-L Representation Approach . As a step toward the goal of analytically 
designing an optimum remote sensing system, the efforts this quarter have 
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been directed towards developing an upper bound on the performance capa- 
bility of such a system to discriminate between two classes. Such per- 
formance bounds will provide guidelines for designing practical systems. 

A preliminary method has been developed to estimate the performance 
bound for some representative field data. The method can be broken down 
into four steps: 1) estimation of the parameters of the process; 2) a 

"whitening" transformation; 3) computation of the likelihood ratio; and 
4) computation of a bound on the performance. 

The estimation of parameters involves finding the maximum likeli- 
hood estimates of the mean and covariance matrices for each class. The 
random process consists of an ensemble of functions of wavelength, and 
the statistics of the process are assumed to be Gaussian. The members of 
this ensemble and their associated probabilities of occurrence are assumed 
to be given for the particular remote sensing problem. 

To simplify the computation of the likelihood ratio and the perfor- 
mance bound, one class is "whitened" with respect to the other class. This 
is done such that the transformed covariance matrix for one class is dia- 
gonal and for the other class the covariance matrix is the identity matrix. 

The optimum classifier is derived from the expression for the likeli- 
hood ratio. The likelihood is defined as the ratio of the probability den- 
sity functions of the two classes [1]. Classification can be performed 
by computing this ratio for each sample function, comparing the computed 
value with a threshold, and deciding in favor of a class depending on 
whether the ratio was greater than or less than the threshold. 

A bound on the performance is obtained by computing the Bhattacharyya 
distance between the classes [1]. The method described above makes the com- 
putation of this distance particularly simple. The bound on the probability 
of error is exponentially related to the Bhattacharyya distance. 

An experimental software system has been set up to implement the four 
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steps of the method and estimate the capability of the system to discrimi- 
nate between any two classes. A trial experiment using wheat and fallow 
data is being used as a design and test mechanism. Approximately thxrty 
sample functions from each of the classes were taken from S191 helicop- 
ter data over North Dakota. Each sample function consists of twenty bands 
over the interval 0.4 to 2.4 micrometers. Use is made of the Karhunen- 
Loeve expansion to generate the transformation vectors for the whitening 
transformation. 

In the coming quarter work will continue to further refine the opti- 
mum system and compare it with practical systems. 

Reference: 

[1] Van Trees, H.L. (1968), Detection, Estimation and Modulation Theory , 
Part I , John Wiley and Sons, New York. 


Information Theory Approach . The approach taken to the spectral 
basis function selection problem is to model the spectral response as a 
portion of a sample function from a stochastic process in wavelength (see 
previous quarterly reports) . It can be shown that between any two spectral 
wavelengths and X 2 the information is given by: 

I(Xi,X2) = % h(X,X)dX 

where h(X,X) is the optimal Weiner filter for estimating the spectral res- 
ponse of the actual scene from the observed spectral process. Evaluation 
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of h(A,A) requires the construction of an adequate model for the spec- 
tral process of a scene. 

It was decided to model the scene in a manner such that Kalman fil- 
tering techniques may be used in the determination of h(A,A) and hence 
the information in a spectral band. 

Several classes of models were hypothesized and several models in 
each class were identified. The identification techniques used were based 
on maximum likelihood identification principles. Of the different classes 
of models considered, the integrated autoregressive models proved to be the 
most promising. 

The best model in the class of integrated autoregressive models was 
found to be: 

Vy(k) = <^iVyk(k-l) + Vy(k-2) + “yVyCk-?) + w(k) 


where 


Vy(k) = yCk) - y(k-l) 

y(k) = spectral response of the scene at wavelength k. 
w(k) = independent, identically distributed Gaussian random 
variables with zero mean and variance p. 

Further w(k) and y(k-j) are assumed independent for all j ^ 1. 

Several tests were conducted to compare the model to empirxcal data. 
These tests were used as a validation procedure for the model. Tests per- 
formed to validate the assumptions above were: 


1. zero mean test 
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2. test for statistical independence 

3. test for absence of sinusoidal components 

Some simulation tests were performed to compare the statistical charac- 
teristics of the model and the empirical data. These tests considered: 

1. correlation similarity 

2. spectral similarity. 

The above model seems to pass the validation tests and hence may be 
considered a model for the spectral process. 

The next step is to put the above model into a form that is amenable 
to Kalman procedures for the determination of h(X,X). Evaluation of the 
information then will follow. It is these tasks that will be addressed in 
the next quarter. 

Task 6: Develop Numerical Methods Classification Error Prediction 


In the previous quarters, an algorithm was developed for estimating 
the overall probability of misclasslf icatlon (PMC) using numerical inte- 
gration technique and requiring the mean and covariance matrices of dif- 
ferent classes. Then the method was applied to the multispectral air- 
craft dater and results compared with the thematic mapper simulation exper- 
iment. These results were reported in the third quarter report. During 
the fourth quarter the scanner modeling effort (Task 7) was pursued and 
further work on Task 6 was delayed. 
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Tasks 8> 9; Implement Multispectral Scanner Model and Evaluate 

The step that has been taken in the quarter is in the direction of 
formulating an input-output relationship for the multispectral scanner 
system. This approach would enable us to analyze the performance of the 
current system; moreover, once put in a proper mathematical/statistical 
model, the parameters of the system can be manipulated to optimize a given 
index of performance. Specifically, it is desired to obtain the PMC for 
a given set of classes once it has passed through the scanner. The com- 
parison of classification accuracy at the input and output of the scanner 
reveals the degradation of the data quality due to additive noise and 
blurring properties of the scanner system. 

Formulation of the problem: 

The scanner is characterized by a 2-dimensional shift- invariant point 
spread function. Independent random additive noise is added at the output. 
See Figure 2.4-1 





m(x,y) 

+ 




Figure 2.4-1. Scanner aperture model. 


f(x,y) is a sample function drawn from a Gaussian, stationary random field 
with separable autocorrelation function. 

Having the statistic of f(x,y) (i.e., mean and covariance matrix) there 
are at least two ways one can obtain the statistic of g(x,y); 1) Matrix 
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operation; 2) spectral analysis. 

1. A convolution in the discrete domain can be carried out 
using matrix operations 


G = H F 

where F is the input matrix, H is the "circulant” PSF 

I 

matrix (i.e., each row is the cyclic shift of the pre- 
vious row) , and G is the output matrix. 

The main drawback here is the sheer size of the matrix 
involved. The mean of G is an NXN matrix and its covar- 
iance contains x .-elements . Even for a modest size 

like N = 10, the covariance matrix will contain 10,000 
elements. Moreover most of the information contained 
in these matrices is not utilized as long as the sta- 
tionarlty and ergodicity assumptions are invoked. For 
example, a scalar would he a valid and sufficient quan- 
tity for the mean of the process under those assumptions. 

2. Co- and Cross-spectral analysis: Here using the estab- 

lished techniques we represent the statistics of the out- 
put process in an algebraic form such that parameter man- 
ipulation is considerably simplified. We attempt to ob- 
tain the autocovariance of the output process, because 
from the standpoint of practical applications, covariance 
and spectral characteristics can be conveniently used 
to get the statistics of the output of a linear shift- 
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invariant system operating in a stationary random environment. Let 


2 2 
X V 

h(x,y) = ^ e"^2a ^ + 2a 

2ira a X y 

X y ■' 


R-(a,3) = e a e a normalize variance = 1 
t X y 


where h(x,y) and R^(a,3) are the scanner PSF and input autocorrelation 
function respectively. Now» it is well known that 


Sg(u,v) = |h(u,v) j^S^(u,v) 
= M(u,v)S^(u,v) 


thus 


Rg(a,3) = R^(a,3)*m(a,g) 

= // R^(a-x,3-y)m(x,y) dx dy 

Carrying out this integration, using the separability of h and R, the output 
sample-to-saraple autocorrelation function will be 

R (a) = e^"“^ q(/~2~ a - M) + Q(/2 a + M) 

y ^ /2 a ^ /T* a 

X X 

where Q(x) is the area under normal density function in the (x,«>) domain. 

R (B) is exactly the same with a, a and a replaced by B, b and a , 

S X y 

respectively. 

It should be mentioned that the entire entry of the aforementioned co- 
variance matrices can be obtained using this expression but new random access 
to its elements is more convenient. In fact, at LARS, the value of interest 
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is the autocovariance at zero lag, i.e., R (0,0). 

Cross-spectral correlation can be similarly computed but the quantity 
required isi the cross correlation at the input. Therefore to obtain the out- 
put statistics, it is essential that a model for the auto-' and cross-correlation 
functions be determined. 


Experimental results: 

Experimental results at this stage, due to data availability, is limited. 
The parameters to be determined are a and b in the autocorrelation function 
expression. The following method is employed: a preferably large field of a 

particular cover type is selected and the mean-lagged product is calculated 
for both horizontal and vertical directions. The maximum lag is chosen to 
be about a quarter of number of samples in each direction. Then a straight 
line is fitted through the natural logarithm of the correlation function and 
the slope is determined using a generalized least-squares criterion. The 
value of this slope is the desired parameter. It should be mentxoned that 
not all cover types exhibit a strictly exponential correlation function and 
deviation from this shape is observed across the fields. 

For Illustration purposes We show an application of the above method 
to a sample field of harvested wheat from the 30-meter simulated thematic 
mapper data. 

Field identification: 

Run No. Cover Type Lines Columns #Points 

75001730 Harvested Wheat 141-153 35-46 156 
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ch 2 

a 

b 

a 

b 

0.54 

0.57 ch 2,6 

0.47 

0.77 

ch 6 

0.54 

0.57 ch 2,8 

0.3 

+1.24 

ch 8 

1 

0.72 

0.38 cb 6,8 

0.45 

+1.00 


Exponential Model Parameters for 
Auto- and Cross-Correlation Functions 



e^ 

I" 

ch 2 

1.82 X 10“^/2.06 X 10~^ ch 2,6 

3.19 X 10“^/2.80 X 10”^ 

ch 6 

1.48 X 10“^/2.18 X 10~^ ch 2,8 

1.08 X 10“^/3.64 X 10"^ 

ch 8 

j 

1.76 X 10“^/4.23 X 10~^ ch 6,8 

2.67 X 10~^/9.06 X lO"^ 


Weighted Least Square Error for 
Horizontal/Vertical Exponential Fitting 


By looking at the plots in Figure 2.4-1 and 2.4-2 and weighted errors, 
it is observed that the correlation functions can be represented as decaying 
exponentials with generally good results, with some exceptions, e.g. , line- 
to-line crosscorrelation between channels 6 and 8. Once the parameters in 
the correlation model are determined as shown above, output statistics for 
any cover type can be computed, with the accompanying classification accuracy. 

The preceding results concerning a particular field of harvested wheat 
is obviously not meant to be a complete and exhaustive model for this cover 
type. This process must be repeated for each crop, encompassing many fields 
and at the end, generating a complete and comprehensive correlation model 
for each channel and across-spectral lands . 

This completes the report of progress on the Scanner Parameter Project 


for the fourth quarter of CY77. 
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Figure 2.4-2. Normalized Crosscorrelation Function versus Lag. 
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2.5 Transfer of Computer Image Analysis Techniques 

^ I 

The successful transfer of computer analysis techniques from a research 
environment to an applications environment by means of a remote computer 
terminal requires a careful blending of hardware and software supported with 
adequate training materials and procedures (Figure 2.5-1). Since its instal- 
lation in late 1972, the JSC/LARS remote terminal has incorporated these 
four components in its remote terminal .operations. As the concept has 
grown and developed the relative emphasis in each of these areas has varied 
according to current needs and requirements. This is shown in Figure 2.5-2. 
This contract year has seen an increasing emphasis on upgrading the hard- 
ware capability of the JSC/LARS terminal as well as the development of 
training materials to test and evaluate ‘new classification algorithms. 

This section summarizes the major activities carried out under Task 
2.5 with emphasis on accomplishments made during the period March 1, 1977 
through May 31, 1977. The major activities were: 

Remote Terminal 

2780 Support 

Data 100 Installation Plan 
LITER Development 

Development of Education and Training Materials 
ECHO Case Study 

I 

Progress in each of these areas are described in the paragraphs below. 
Following this is a proposed technology transfer milestone chart and a 
discussion of work planned for the future, 

2780 Support . Support of the JSC/LARS 2780 remote terminal Includes 
maintaining at the LARS computer facility the necessary hardware and soft- 
ware, providing computer services, providing hardware/software system consul- 
tation, and developing concepts for an improved remote terminal system. 



REMOTE TERMINAL SYSTEM 



HARDWARE 
SOFTWARE 
TRAINING MATERIALS 
PROCEDURES 


Figure 2.5-1 
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Year 

1972 

1973 

1974 

1975 

1976 

1977 


Figure 


Major Activities 


Emphasis 


Installation of remote terminal 


LARSYS 3.0 release. Revised 
training package 

LACIE Analyst training 


LARSYS 3.1 release, LANDSAT 
Case Study 


Simulation exercises, train- 
ing procedures documented 


Terminal upgrade, ECHO Case 
Study 


Hardware, training, procedures 
Software, training 

Training 

Software, training 

Training, procedures 

Hardware, new algorithm train- 
ing 


2.5-2 Historical picture of JSC/LARS remote terminal activities 

showing how relative emphasis shifts between hardware, soft- 
ware, training and procedures depending on current needs. 
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Hardware support consists of maintaining communication modems at LARS 
and a port on the 3705 communication system. Software support includes 
programming the 3705 to recognize the JSC terminal, CP-67 and LARSYS soft- 
ware, and maintaining a virtual machine dedicated to the JSC terminal. 

During this period, communications with remote terminal users at JSC 
increased markedly. It was decided that certain routine requests, such as 
ID changes and data reformatting requests could go directly through Bob 
Peery or Jim Tyler at JSC, with copies of information sent to John Sargent. 

A current list of computer IDs and passwords was mailed to Jim Tyler, and 
several computer IDs have received increased disk storage space during this 
quarter. In April, a memorandum discussing the scheduled Labor Day instal- 
lation of an IBM 370/148 computer here at IjARS was sent to all remote terminal 
sites. The effect of this change on remote terminal users, aside from a few 
days of down time, should be minimal. No hardware changes are required and 
user Interaction with the system will remain almost the same as now. In 
addition, the remote sites will gain the capability to manipulate their 
printer output files. 

At the request of a remote terminal user at JSC, PL/I was made avail- 
able to JSC and has been successfully used. Also during this quarter a LARS 
update to the CP-67/CMS manual was mailed to all remote sites, and a copy of 
additional system information was sent to JSC at the request of a JSC terminal 
user. Two data runs were entered into the LARS data library during this 
period, although one data entry was delayed because one of the CCT’s received 

I 

by LARS was unreadable. 

We have noticed that one or two JSC users are making frequent use of 
the dial-up lines into the Purdue/LARS computer. If this dial-up access is 
desired by JSC for the future, we recommend investigation of obtaining dial- 
up access to the Codex multiplexer at JSC which would then make use of the 
existing leased line between JSC and Purdue/LARS. 
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Computer usage for the period February, 1977 through April, 1977 was 

i 

as follows; 


Month 

CPU hours used 

Hrs. terminals attached 

February 

12.2 hours 

291 hours 

Ifarch 

17.1 hours 

378 hours 

April 

11.4 hours 

289 hours 


Data 100 Installation Plan . In response to NASA’s request (December 
20, 1976), lAES has estimated the performance, cost, and schedule for instal- 

t 

la t ion of a tape drive capability on the JSC/LARS terminal. During the 
December 1, 1976 to February 28, 1977 quarter three approaches were inves- 
tigated. Their throughput capabilities, hardware cost and requirements, 
software requirements and leadtimes were presented in the December 1, 1976 - 

i 

February 28, 1977 quarterly report. Subsequent to the writing of that report 
a fourth option utilizing a Data 100 terminal became feasible. This possi- 
bility was discussed with JSC personnel during the March 2-3, 1977 quarterly 
review. The positive response to this option resulted in the preparation of 
a plan for installing a Data 100 remote terminal at JSC. This plan, which 
is presented below, was transmitted to Mr. Don Hay on April 1, 1977. 

o o .0 0 o o 

PLAN FOR INSTALLATION OF DATA 100 REMOTE TERMINAL 

Howard L. Grams 
Terry L. Phillips 

Laboratory for Applications of Remote Sensing 

Purdue University 
March 1977: 

This document is intended to present specific details of the 
proposed implementation of hardware enhancements to the existing 
remote terminal located at the Earth Observation Division at the 
Johnson Space Center. The hardware currently includes an IBM 
2780 system and is connected to the Eurdue/LARS computer. The 
terminal can be upgraded to allow functions to be performed that 
now cannot be performed at all, and to allow functions that can 
be performed now to be performed faster or more conveniently. 
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The proposed hardware is expected to be significantly more 
reliable than the present hardware. The "human factors" aspect 
of using the new hardware should be improved from that of the present 
hardware, and the retraining aspect is expected to be minimal 
to perform present functions. 

The proposed change is to replace the present IBM 2780 terminal . 
with an upward compatible terminal. At present this is expected 
to be a Data 100 Model 78 as described in the subsection entitled 
"Tentative Hardware Budget". The new terminal would have a magnetic 
tape capability that could be used either as a source for data 
sent to the lARS computer or a sink for data received from the 
LARS system. The new terminal would otherwise be compatible with 
the existing 2780 (it can read cards, print output, and punch cards). 

In some cases these functions would be faster than those of the 
2780 and also more convenient frorri a, human factors point of view. 

In addition, the existing Codex 7200 modems and 800 multi- 
plexors would be replaced with newer hardware that allows higher 
speed transmission. Hov;ever, the existing C2-conditioned phone 
line would still be used. 

The benefits and purposes of the new terminal are included 
in this report as well as a milestone plan for installation. 

Benefits and Purposes 

The proposed implementation is an amalgam of need, ideas, 
wishes, and suggestions fiom both JSC and LARS personnel. Its 
benefits include (but are not limited to) the following: 

A. Allow JSC personnel to send test data to LARS so they 
can have it immediately available to use for research, 
test, and evaluation efforts using the LARS computer 
hardware/software systems. ^ 

B. Provide means for JSC personnel to test and evaluate 
the ECHO technology presently under development at LARS. 

This would be facilitated by enhanced access to the 
algorithm implementation on the LARS computer as well 

as by the capability to quickly and conveniently send 
test data. 

C. Improved access to the Field Measurements data base 
for EOD personnel at JSC. 

D. Enhance the present operational procedures at JSC for 
CAMS/CAS by providing a more convenient medium (magnetic 
tape) for the operation than presently used (punched 
cards) for data read into software in the LARS computer 
and output from that software. 
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MILESTONES 

This discussion is centered around a possible proposed 
scenario, with certain assumptions made concerning starting dates, 
necessary decision dates and so forth. A milestone chart for this 
scenario is attached. It should be noted that any delay of each 
Particular milestone will produce a corresponding length delay 
in all other milestones downline. 

DATE 

RESPONSIBLE 

DESCRIPTIONS 

May 1 

Purdue 

Submit proposed hardware specifications 
to JSC. This includes physical site re- 
quirements information. 

June 1 

JSC 

Make final decision to proceed with 
project. Approve funding arrangements. 
Approve the choice of hardware recommended 
by Purdue for procurement. 

October 1 

JSC 

Order hardware for the terminal and 
modems by July 1. Since delivery time is 
expected to be 90 days, any delay will cause 
the final installation date to slip a corres- 
ponding amount. 

Order, power installation and physical 
site preparation. 

Install power receptacles and have 
physical site for new terminal hardware 
ready by September 15. 

Install new terminal hardware, modems, 
and related equipment on October 1. 

August 1 

JSC 

Approve plans presented on July 1 
by Purdue : for software development to 
support the new tape drive capability. 

August 1 

JSC 

Approve plans presented on July 1 
by Purdue for retraining terminal users. 

November 1 

Purdue/JSC 

Complete retraining of users so that 
2780 can be removed. Complete testing of 
new hardware to insure it performs all 
current functions. Complete testing and 
declare operational the new functional 
capabilities that utilize the magnetic 
tape uniti' Train users at JSC in its use. 
Remove 2780 hardware. 

November 15 

Purdue/JSC 

1 

Meet to evaluate the results of the 
project and the state of implementation 
of the functions using the new magnetic 
tape hardware. Identify problem areas 
for future effort. 
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Tentative Hardware Budget 

The proposed new hardware to replace the existing IBM 2780 
terminal is a Data-100 Model 78 terminal with a magnetic tape 
capability. Although this hardware is available for lease (12 
months minimum, month- to-month thereafter) or for purcahse, 
it is recommended to lease since LARS expects to develop even 
more cost effective hardware in the reasonably near future. 


The hardware suggested includes; 


78-111 

Terminal Control Unit 

$ 355/mo 

8801 

IBM 2780 Emulation Feature 

26/mo 

8130 

Card reader adapter 

13/mo 

8132 

Card punch adapter 

102/mo 

8140 

Line printer adapter 

40/mo 

8190 

Magnetic tape adapter 

82/m6 

78-202 

9- track 1800 bpi Magnetic tape 

219/mo 

8878 

4K MT Blocking feature 

32/mo 

78-304 

150 CPM Card Reader 

109/mo 

78-302 

Card Punch 

288/mo 

78-405 

400 LPM Line Printer 

596/mo 

8442 

132 Column Option 

77/mo 

78-702 

RS-232 Synchronous Communications 

98/mo 

8721 

9600 BPS option 

40/mo 


TOTAL $2, 077/mo 

The proposed modem configuration consists of; 

2 Codex LSI-96 Modems @$185/mo $370/mo 

1 Codex 910 multiplexor with RAM 124/mo 

1 Codex 910 multiplexor with dual RDM 114/mo 

Miscellaneous Options (Mounting HW, ' --20/rao 

cabling, etc.) 


$628/mo 

Dual Configuration RDM for multiplexor $180 

. (purchase) 
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We are currently waiting to receive a response from JSC with regard to 
this plan. 


LITER Development . The concept of an improved remote terminal based 
around a minicomputer supported by an array of special purpose and general 
purpose peripheral equipment was described in the June 1 - August 31, 1976 
quarterly progress report. Progress on the development of this terminal, 
known as the LARS Intelligent Terminal for Earth Resources (LITER) , during 
the past year has been limited. This is due to the relatively low funding 
level and, as reported in the September 1 - November 30, 1976 quarterly 
report, a decision, based on discussions ' with JSC, to lower its priority in 
favor of raising the priority of proposing solutions to JSC's expressed need 

I 

to have an improved capability for transmitting data over the remote terminal 

I 

link . ‘ 


Progress which has been made during the past year falls into three cate- 
gories — planning, equipment installation, and capability demonstrations, 

f 

Planning — Implementation of the concept has been planned 
in three phases: I - a developmental prototype, II - imple- 

mentation of user oriented software, and III - hardware 
additions and further enhancement of the system software. 

Only a portion of Phase I has been completed this year. 

Equipment Installation — The following equipment asso- 
ciated with the LITER system has been installed and 
tested: 

PDF - 11/34 Minicomputer 
Varian 4211 Electrostatic Printer/Plotter 
Tales 660B Table Digitizer 
Info ton Vistar/GTX CRT terminal 
plus necessary interface equipment. 

Capability Demonstrations — Progress in demonstrating 
the capabilities of the LITER system has been limited. 

Each of the major capabilities of the system have been 
tested and sample products produced but little progress 
has been made in documenting the system capabilities. 
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Development of Education and Training Materials . During the past year 
there has been a shift in emphasis with regard to the development of education 
and training materials under SR & T support. Effort on materials dealing 
with fundamental concepts and current application tasks has been decreasing 
and an increase in effort has been made in the development of materials to 
support the transfer of new technology to JSC. This reorientation of prior- 
ities is consistent with the proposed plans to upgrade the JSC/LARS remote 
terminal to make it a more effective tool for the test and evaluation of new 
analysis procedures. 

A) Items completed. LARS Information Note 052977 entitled "The FOCUS 

Series: A Collection of Single-Concept Remote Sensing Educational Materials" 

was published this past quarter. This Information Note is a compilation of 
all titles in the FOCUS series, each of which consists of a diagram or photo- 
graph and an extended caption of three to four hundred words treating a 
single concept. Purposely intended to be brief, individual titles in the 
series do not lend themselves to being submitted as technical reports. In 
order that the existance and rationale behind the series may be disseminated^ 
the titles currently available have been' collected in this Information Note. 

B) Items in progress. An analysis case study is being developed in 
order to provide education and training in the use of the supervised ECHO 
classification algorithm. The case study is scheduled to be completed with- 
in three months after the release of the ECHO software. Unlike the case 
studies developed for 2780 remote terminal training the ECHO case study is 
being written so that it is not dependent upon a particular hardware conflg- 

r 

uration. This will give the case study wider utility and will facilitate 
the test and evaluation of the ECHO processor as it might be implemented on 
JSC’s computer system. The new format essentially simulates a batch mode 
computer environment. The case study is divided into a number of steps 
corresponding to the analysis sequence. The student will interact with a 
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tutor at the conclusion of activities stratigically placed through the 
materials. Appropriate computer printouts will be available at key points 
simulating a batch machine operation. 

Two data sets have been selected for use in the case study. One will 
be used to provide a detailed analysis example, the other will be used in a 
series of exercises. Persons working the exercises will independently com- 
plete an ECHO classification using available computer facilities. 

Each step in the analysis will be discussed at three levels; an overview 
describing the purpose, objectives and techniques used; a second level dealing 
with the specifics of a sample analysis; and, a third level which will allow 
the student to carry out his own ECHO analysis. A draft of the case study has 
been completed and is presently undergoing internal review. 

Proposed Technology Transfer Milestones . Figure 2.5-3 shows a proposed 
set of technology transfer milestones through 1980. The major milestones 
include installing and supporting a Data 100 terminal at JSC, and demon- 
strating the feasibility. Installing and supporting a LITER terminal at JSC. 
These long term objectives are all aimed at providing a significantly 
improved remote terminal which can serve the research, test, and evaluation 
needs of JSC. 

Plans for Next Quarter . Plans for activities under this task during 
the next quarter center around the Data 100 Installation Plan milestones 
(see page 2.5-5 through page 2.5-9). In addition to the items specifically 
mentioned there the ECHO Case Study is scheduled to be completed during the 
next quarter. Its use will be incorporated into the Data 100 terminal train- 
ing plan in order to demonstrate how the remote terminal system can serve 
as an RT & E tool. Proposals for using the terminal to evaluate Procedure 1 
(PI) will also be part of the Data 100 training plan. 
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Problems . Purdue/LARS has not received any response from NASA/JSC con- 
cerning the Data 100 Installation Plan. Delay beyond June 1, 1977 will 
cause a corresponding delay in installing the terminal. 



Figure 2.5-3 
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2.6 Large Area Crop Inventory Design 

In March the Large Area Crop Inventory Design (LACID) task at the 
request of NASA/ JSC began to concentrate Its activities on the require- 
ments for a food and fiber information system. On March 30-31, R.B. 
MacDonald and three of his staff met with staff and several department 

r 

heads of the Purdue School of Agriculture to explore possible roles of 
the School and LARS in assisting JSC to define the information require- 
ments of a global food and fiber information system. 

As a result of these and later meetings involving the School of 
Agriculture administration, the Implementation Committee was encouraged 

I 

to proceed with organizational and work plans for implementation of the 
study. 

Four staff members of the School of Agriculture are meeting regu- 
larly to implement the study. This group,' known as the Implementation 
Committee, is composed of Prof. Marion Baumgardner, Chairman, Agronomy; 
Dr. Marvin Bauer, Agronomy /LARS ; Marshall Martin, Agricultural Economics; 
and Robert Peart, Agricultural Engineering. An advisory committee of 
20 scientists from nine departments has been named. This committee will 
provide counsel for Implementation Committee and review materials pre- 
pared by the Implementation Committee. Several potential consultants 
with experience in agricultural policy, rangeland and forestry have been 
contacted. 

At the present time we are acquiring 'and reviewing reference mater- 
ials and reports from previous studies which may have a bearing on the 
definition of a global food and fiber information system. The proposed 
schedule of activities calls for completion of detailed proposal of 
objectives and approach by mid- June. Plans are being made for a work- 
shop in September with participants from the scientific community, agri- 
cultural industries, and international development organizations. 
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2.7 Forestry Applications Project 
Introduction 

This task activity has been directed toward the development of inventory 
schemes which utilize Landsat data and computer-aided analysis techniques . 

The ultimate objective is aimed at developing various approaches for using 
Landsat data to conduct inventories of the 1.6 billion acres of forest and 
range land resources in the United States. Emphasis is being placed on 
effective, statistically sound approaches which a-priori inputs, such as 
climate, topography and previous inventory information. 

Background 

The approach has developed through several phases, each becoming better 
defined and structured than the phase proceeding. The initial phase or 
activity dealt with a review of forest inventory literature which emphasized 
sampling techniques developed since the launch of Landsat 1. The best focus 
for the project was acquired through a review of the 1975 assessment of the 
Resources Planning Act completed by the U.S. Forest Service. The result of 
the review is summarized in Table 2.7-1. Specific needs were identified that 
could be addressed through the application of remote sensing technology. 

During the second phase our attention was focused on the forestry sector 
of the Resources Planning Act. Specifically, we worked to define how a 
Landsat analysis system could benefit the ungoing forest survey activity 
conducted by the Forest Experiment Stations. In Effect the task objective 
had been narrowed from broad-based large-scale surveys to regional single 
purpose surveys. By so doing, we were able to define benchmarks (forest 
survey standards) against which we could assess the effect of Landsat inputs. 
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Table 2.7-1. Summary Review of Resource Planning Act Assessment (1) . 


System 


Unit of Measurable Beneficient 

Measure by R. S. to Decision 

(2) (3) (4) 


Timber 


Type 

Area 

+ 

Ownership 

Area 

0 

Productivity 

Volume 

0 

Range 

Type 

Area 

+ 

Ownership 

Area 

0 

Productivity 

Volume 

0 

Water 

Typ e-Impound ed 

Area 

+ 

Flowing 

Area 

0 

Use 

Volume 

- 

Fish and Wildlife 

Habitat Type 

Area 

+ 

Population Estimates 

Census 

- 

Recreation and Wilderness 

Participation 

Census 

- 

Human and 

Community Development 

Economic 



D 

I 

I 

D 

I 

I 

D 

I 

I 

I 

I 

I 

N 


(1) ref: The Nation’s Renewable Resources-An Assessment, 1975, USDA- 

Forest Service. 

(2) General classification of units of measure reported in (1) . 

(3) Measurable by remote sensing techniques. 

Includes all measures possible from satellite and/or aircraft 
systems . 

+ Directly measurable from satellite or aircraft data alone. 

0 Measures can be inferred with the support of ancillary type 
data. 

- Measure which is not even indirectly affected by satellite 
or aircraft inputs. 

(4) Benefit of Remote Sensing inputs to Decision Process . 

Inferred based on Inputs . 

D - Direct not utilizing other Inputs, 

1 - Indirect utilizing other inputs. 

N ~ Not Obvious. 
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Forest Survey is required by congress to meet certain pre-defined stan- 
dards for measuring the areal extent (expressed in acres) and quantity 
(expressed in cubic feet) of the available resource. Our study indicated that 
most sampling schemes considered quantity of the resource a critical variable, 
while area, which is the most easily measured parameter from Landsat, is all 
but ignored. Therefore, most of these approaches require complex statistical 
manipulation of various data sources to identify field plots which are necessary 
in order to acquire volume. 

Forest Survey, depends on estimates of forest acreage in order to allocate 
field plots. If we assume that a Landsat system can accurately measure areal 
extent of forest, then we may anticipate some improvement over the current 
approach. The second phase was directed toward formulating assumptions which 
could be used to answer the question relative to Landsat 's ability to provide 
the acreage inputs necessary for Forest Survey. 

In addition, we are preparing to Investigate the utility of the Forest 
Survey permanent remeasurement plot network to the design. Plot data is being 
collected for sites in Indiana where Landsat data is available. The relation- 
ship of field measured parameters to spectral response will be studied. Finally, 
an assessment of the utility of this data to the design problem will be made. 

Phase three of this task will identify recommended approaches to the inven- 
tory problem. 

Approach 

After review of the Resources Planning Act Assessment of 1975 and Forest 
Survey requirements, we outlined a general approach for the use of Landsat 
data. The approach is graphically displayed in Figure 2.7-1. Table 2.7-2 
details the anticipated use, potential benefits and detail of ancillary input 

I 

required for each level of classification intensity. 
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Figure 2.7-1. Schematic flow anticipated in design definition. 
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Table 2.7-2. Anticipated use, potential benefits, and ancillary i 
classification intensity. 


Design 

Intensity 


Input from 
Landsat 
Analysis 


Unit 


Level of Use 


1st level Forest 


Non-Forest 


Acres - Proportional allocation 
of points for intensive 
airphoto and ground 


2nd level Non-Forest 

(< 10% cover) Acres - Same as level 1 


Commercial 

Non-Commercial 


- could optimize plot 
allocation based on 
previous inventory 


3rd level 


Information 

Partitions 


Acres 


- Same as level 2 


Within Commer- 
cial Forest 


inputs for the various levels of 


Potential Benefits 


- Timeliness 

- Speed of Analysis 

- More Objective 


- Same" as -level 1 

- Minimize temporary 
plot location and 
field work 

- Improve allocation 
of field work 

- Increase efficiency 
of survey 

- Decrease frequency 
between surveys 

- Yearly monitoring of 
forest change 

- Same as level 2 

- Improve efficiency 

- Improve yearly pre- 
diction of available 
resource supply 

- Monitor forest change 
in greater detail 


Ancillary Data Needs 


- Maximum ancillary 
data needs (air 
photo & ground) 

- Minimum input, only 
acreage 

- Maximum photo and 
field plots 


- Possible decrease 
in field plots 

- Probable reduction 
of photo plots 


- Minimize amount of 
ancillary Inputs 
(fewer photo points, 
no temporary plots) 
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This approach assiimes that Landsat data can be used to delineate areal 
extent of forest test areas, consisting of counties within the Great Lakes 
Watershed will be used as benchmarks against which Landsat classifications 
will be compared. Data are currently available for Michigan and Wisconsin, 
and are being acquired for Minnesota and New York. Statistical tests will be 
run on the Landsat and Forest Survey results. 

The efficiency, or timeliness, of obtaining classification results from 
Landsat are of concern for the Inventory design. Any advantage of a Landsat 

•t 

based inventory system may be greatly diminished if the time required to pro- 
duce the acreage estimates was greater than or equal to current procedures. 

A small study has been defined to provide quantative information on 
classification time. Carlton County, MN i’s being classified using both a 
maximum likelihood approach and a levels slicing approach. The "levels" 
approach offers a potential savings in classification time given the analyst 
has sufficient or prior knowledge to set the levels. The time required to 
achieve similar classification accuracies for a fores t/non-forest breakdown 
will be composed. Recommendations regarding the "levels classify" approach as 
part of a multi-level sampling strategy will be made in the final report. 

The second level of classification intensity is intended to optimize a 
prior knowledge about the forest resource through use of Forest Survey per- 
manent sample plots. Figure 2.7-2 shows a portion of a field plot sheet used 
to record permanent plot data. Approximately 30 such plots are available for 
two sites in Indiana. These data are currently being studied to determine 
how they can best be used. 

The last level of classification intensity, the level of maximum infor- 
mation detail, will not be addressed by this study. Level 3 knowledge assumes 
that a high correlation exists between spectral categories and useable infor- 
mation. The investigation into this level of detail is beyond the scope of 
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this task. Based on our understanding of the scene components and the current 
technology we will be prepared to make recommendations regarding future areas 
of study. 
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Figure 2.7-2 Portion of forest inventory field record showxng the type of information collected for 
pemnanent field plots. 
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2.8 Regional Applications Project 
Objective 

The overall objective of the Regional Applications Project was the 
development of a computer procedure for (1) Inventorying and categorizing 
Texas coastal zone environments and (2) monitoring and identifying signifi- 
cant alterations in these environments. 

Within this objective was the goal of identifying and delineating 
primary level coastal zone features based upon the association between 
Landsat multispectral data and known ground truth information. On site 
inspection, interpretation of aerial photographs, and Texas Bureau of 
Economics Geology maps were to form the ground truth information. 

Approach 

Previous research by Purdue/LARS involved the classification of Landsat 
MSS data to produce an inventory of environmental units within the Texas 
Coastal zone. For the 1976/77 contract year, the emphasis was to be on the 
development, refinement, and documentation of a change detection procedure. 
This effort was to incorporate the comparison and analysis of Landsat MSS 
data accumulated over different time periods for the Matagorda Bay Coastal 
Zone study site (three USGS 7% minute quadrangles, i.e.. Pass Cavallo, Aust- 
well, and Port O'Connor), In the development and refinement of a change 
detection procedure, the following efforts were to be accomplished by Purdue/ 
LARS: (1) Determine the extent to which detectable seasonal changes in the 

coastal zone classifications could be identified, (2) Determine the extent 
to which detectable permanent changes in the coastal zone classifications 
could be identified. (3) Determine methods for delineating between seasonal 
and permanent variations for specific coastal zone features. (4) Register 



2 . 8-2 


detected change classes to 1:24,000 USGS quadrangles using existing LARS 
capabilities. (5) Define significant changes to be established by the 
user agency, i.e,, (State of Texas). Using these definitions a measure of 
the degree of accuracy of the change detection procedure was to be made. 

(6) Transfer the change detection procedure with software, algorithms 
and user documentation to NASA/ JSC. 

Accomplishments 

The accomplishments under this task, a six-month investigation ter- 
minated by the sponsor on 30 November 1976, are summarized as follows. 

Of the four change detection techniques investigated during the 
past contract year (post classification comparison, delta data, 
spectral/ temporal, and layered spectral temporal), the post classification 

I 

comparison technique was selected fcjr further development. This was -based 
upon test performances of the four change detection methods, straightforward 
ness of the procedures and the output products desired. 

A standardized "modified" unsupervised classification procedure for 
analyzing the Texas Coastal Zone data was compiled. This procedure was 
developed in order that all quadrangles in the study area would be classi- 
fied using similar analysis techniques to allow for meaningful comparisons 
and evaluations of the classifications. Single date classifications were 
completed for the three quadrangle areas for the February 25, 1975 date. 
These classifications were to be used as the base data to which classifi- 
cations developed from the later three dates were to be compared in order 
to analyze the post classification change detection technique and to develop 
a standardized change detection procedure. 

A "modified" unsupervised approach was used to classify the Austwell 
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quadrangle. The Spectral Environmental Overlay in conjunction with February 
1975 CIR photography was used to choose training areas. The resulting clas- 
sification contained 35 spectrally separable classes which appeared to char- 
acterize the scene completely. 

Since the Spectral Environmental Overlay contained only 25 classes, a 
one-to-one correspondence between the overlay classes and the 35 spectral 
classes could not be established. Preliminary comparisons show that some 
classes of the hierarchical classification are represented by more than one 
spectral class, and there are some hierarchical classes (e.g,, TFi, WB) that 
are not represented in the spectral classification of the Austwell area. 

The absence of these classes can be explained either by there being so few 
data points available to use as training samples to represent the classes 
or the classes are not spectrally separable from other classes. It appears 
that the latter case may prevail. Grouping of certain members of the 35 
spectral classes was done to gain the maximum correlation between the spec- 
tral classes and the hierarchy classification. 

Thus, an additional task involved the modification of the spectral 
hierarchy developed by JSC/LEC. The original hierarchy was developed using 
photo interpretive information such as texture, spatial distribution, sub- 
jective cultural Interpretations and spectral response. Since the Landsat 
multispectral scanner records data predominantly spectral in nature, any 
hierarchy to be used with scanner data must be spectrally derived. Table 1 
lists the 35 spectral classes Identified in the Austwell quadrangle area 
and their corresponding JSC/LEC hierarchical classes. The final version 
of the modified spectral hierarchy was not completed. 

Another area of involvement consisted of investigation of methods for 
improving the geometric accuracy of the registered coastal zone quadrangles. 
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In the Port O'Connor quadrangle there are no control points in the lower 
right quarter since this is an offshore region, A one-to-two line geometric 
stretch error is observed and this is thought to be due to the effect of 
inaccurate quadratic terms caused by poor point distribution. The original 
geometrically corrected Port O'Connor quadrangle data were inspected to 
obtain all possible ground control point locations. The same points were 
located in the uncorrected data and on the USGS map used as a geometric 
reference. The mathematical functions tested were the affine and the col- 
linearity geometric corrections . The results of these methods were compared 
to the least squares bi-quadratic method used in the original corrections. 
The methods were tested by selecting check points in the lower water covered 
portion of the Port O'Connor quadrangle where no control points existed. 
Results indicated, that the affine geometric correction function is the most 
accurate of the two tested in areas where control points are not well dis- 
tributed over the area to be corrected. 

The RAP simulation study was completed during this period. This study 
consists of a documented procedures package for classifying the Pass Cavallo 
quadrangle in the Texas Coastal Zone. This package will be valuable in 
transferring analysis technique procedures, as applied to coastal zone fea- 
tures, to persons new to data analysis or coastal zone environments. 

Technical Problems Encountered 

The most serious difficulty encountered during this effort was the 
lack of available data for use in developing procedures. 

Recommenda tions 

Work should continue on the development of procedures for Inventorying 
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and monitoring the Texas Coastal Zone. However, for this effort to be 
successful, a positive effort must be made to collect adequate reference 
data coincident with the desired Landsat overpasses . Without coincident 
reference data performance criteria cannot be properly evaluated. 

Even though the post classification comparison change detection tech- 
nique was selected for this portion of the study, it is felt that the 
layered spectral/temporal change detection procedure holds great promise 
for monitoring change. This technique should be reevaluated when an appro- 
priate data set (coincident reference and Landsat data) has been assembled. 
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SPECTRAL CLASS DEFINITION - AUSTWELL QUADRANGLE 


Spectral 

Class 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 


Class Definition 


Hierarchical Class (es) 
Represented 


17, 18 
19, 20 

21 


22 


low vegetative cover on light (bright) soil He, Heg 

sparse low vegetation Het 

mixed woody /herbaceous cover-f airly green Mdg , HGg, some Md. 

wet herbaceous-lush cover He , some Het, Heg. 

herbaceous -covered with water Hei 

herbaceous-complete inundated Hei 

mixed-woody and herbaceous-f airly dense He, Wf 

bright bare soil-wet spread all over in the Bst 

delta 

herbaceous-wet-occurring mainly in the delta Hei 

standing water,, dark-low reflectance-some some Hei 

vegetation 

herbaceous-farily dense He 

herbaceous emergent wet Hei, Heg 

sparse wet herbaceous-low spectral response Hei, He 

sparse wet herbaceous-low spectral response- , Hei, He 
few points occurring alongside class 13 

bare soil wet Bst 

water in pond next to the refinery (?)- MMw 

high reflectance in class 2 Indicating high 
turbidity-or some pollutant-or red suspended 
sediment 

Sea water classes W 


bright soil (sandy slightly inundated-yf ringing Bs, Ws 

around shallow water ponds and in the delta 

sparse vegetation on light colored soil that is Heg 

wet-few points in the delta 

mixed pixel response-bare soil and wet vegetation 

scattered in the lower delta region and along the 
coast 


23, 24 
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Spectral Hierarchical Class (es) 

Class Class Definition Represented 



25 

"edge" ’Mater class-occurring along the coast 
and fringing shallow water areas in the delta 

Ws 

29, 

30 

bare soil tilled and wet 

Bst 


31 

bright wet surface (concrete) delineating 
the refinery 

MMu 

32, 

33 

very few points picking up the commercial 
in the city residential class 

MMu 

34, 

35 

sea water in the bay area 

W 

26, 

27 

shallow water 

Ws 


28 
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2.9 Interpretation of Thermal Band Data 
Revision of Implementation Plan 

This is a continuing project in a contract extension and a 
quarterly report will be given here in lieu of a previously planned 
final report. The implementation plan is to be revised to this 
extent with subsequent Impact on the milestone chart. 

Major Activities . A simplified radiation model intended to util- 
ize the experimental data acquired previously has been studied in 
this quarter. The objective was to produce a model that would pre- 
dict the effective radiant temperature of the canopy under study, in 
this case North Dakota Spring and Winter Wheat, as a function of 
view angle (8 where 0’ and <j)' are the viewing zenith and 
azmuith angles respectively. The radiant temperature would presumably 
be a function of canopy geometry, solar input power, and soil type. 
With proper agronomic characterization the canopy geometry can 
ultimately be correlated to the yield of the crop. 

Generally two basic approaches could be used: a rather complex 

model that relies on easily produced experimental data or a simple 
model that depends on a relatively complex experiment. The latter 
was chosen because of the difficulties encountered by many investi- 
gators in mathematically describing a complex structure like wheat. 

The model chosen uses a one-layer structure with the heated soil 
at the bottom and the tips of the plants at the top. The quantities 
needed to support the model are the radiant temperature profile T(z), 
where z is the vertical distance in the canopy (soil surface at z = 0) 
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and the view factor P(z,6’,^*) which describes the probability of 
viewing the sky at any height z in the canopy in the direction 
The view factor is directly relatable to the canopy geometry and 
indirectly to the canopy yield. If continuous functions were available 
the radiant temperature of the canopy in direction (0',<j)’) would be 
given by 

z 

f,(T (9', <!>’)) - / P(z,0',<f.') f (T(z)) dz. 

^ V Q J- 

fj^ is derived by applying radiation heat transfer principles to T(z) 
and Ty(0’,(|)') is derived by inverting f 2 which is also based on 
radiation heat transfer principles. In practice, however, continuous 
functions are not available and the model uses discrete approxima- 
tions to the above functions. 

T(z) is obtained using procedures described in previous quarterly 
reports. The experimental setup is schematically described in 
Figure 2.9-1 and is shown photographically in Figure 2.9-2. Some 
typical data products are shown in Figure 2.9-3. Note that the 
thermal scanner is arranged to scan vertically through the canopy and 
the line scan data shown in Figure 2.9-3 B and D, when calibrated, 
directly yields T(z). In all cases the radiant existance is interpre- 
tated as emitting from a target with an emittance of unity prior to 
the calculation of the radiant temperatures. 

The view factor P(z,0 ’,({)’) is obtained experimentally using pro- 
cedures described in LARS Information Note if 120776 CA Laser Technique 
for Characterizing the Geometry of Plant Canopies; Vanderbilt, Silva 
and Bauer) . The procedure Involves the utilization of a laser to 



Down- row View 



Grass 


Figure 2.9-1 Arrangement of Thermal 

Scanner and other instru- 
mentation for Plot #124 
Observation. 




Cross-row View 
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(A) Thermistor Stake in Wheat Canopy 



(B) Arrangement of Instrumentation 


Figure 2.9-2 Photograph of the Thermistor Stake and Instrumentation Arrangement 

(Plot #124, 17 June 1976) 
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(A) Thermal Image, (B) Line Scan, 

full sunlight full sunlight 

(local time 1425) 




(C) Thermal Image, 

thick cloud cover 
(local time 1428) 


(D) Line Scan, 

thick cloud cover 



Figure 2.9-3 Side Thermal Images and Line Scans 
(Plot //124, 17 June 1976) 
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(A) Canopy Side View 



(B) Digitized Thermal Image 


Figure 2.9-4 Photograph and Digitized Thermal Image, Side View of Canopy 
(Small Grain Test Plot #7, 10 July 1976, local time 1801) 
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optically probe the canopy. The laser is aimed into the canopy 
at several zenith angles C®*)* Initially azmuith (<)>’) effects have 
not been considered. In order to adequately characterize the canopy 
sufficient laser data must be acquired. If data are taken over a 
range of zenith angles from 5° to 75° in 10° steps eight groups of laser 
data will result. For at least 1% quantization accuracy 100 laser 
hits are required in each angular group of data. To account for 
canopy variability at least three replications of the data are needed 
which means that 8 x 100 x 3 = 2400 laser hits are required. The 
laser technique is intended to be utilized with an automated system. 
Budget restrictions forced the use of a manual system and, unfortunately, 
insufficient data were take to adequately characterize the canopy. 

As a result the radiation model cannot be fully evaluated with 
available data sets. 

An experiment has been planned for the summer of 1977 in conjunc- 
tion with the Field Measurements Research project that uses a modi- 
fied laser technique, basically manual, that should yield sufficient 
data to characterize the canopy. A backup photographic technique 
will also be used in case the proposed laser technique proves to be 
too cumbersome. The thermal measurements will be repeated in one 
series of experiments during July, 1977, at the Wllliston, N.D. test 
site. 


Technical Problems and Solutions . Insufficient laser data were 
taken to adequately characterize the wheat canopy under study. A 
revised procedure has been developed for a limited series of 
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experiments in July, 1977, that should yield sufficient data to 
support the radiation model. 

Plans for the Next Quarter . Additional experimental data will 
be taken to support the radiation model. 
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2,10 Super Site Data Management 
Major Activities 

The quarter's activities include continued work towards verifying 
and documenting the quality of the truck-mounted and helicopter-moiinted 
spectrometer data. A Field Measurements Spectral Data Verification 
Report is being written which, discusses the quality and the correlation 
of the data for the major sensors. The sections planned at present for 
the report include 

1) Spectrometer Calibration Method 

2) NASA/JSC Field Spectrometer System CFSS) Data 

3) NASA/JSC Field Signature Acquisition System (FSAS) Data 

4) Purdue/LARS Exotech 20C Data 

5) NASA/ERL Exotech 20D Data 

6) Spectrometer Systems Data Correlation 

7) Aircraft Multispectiral Scanner Data 

FSS(S-191H) Spectrometer Data 

The primary verification '.work for the FSS data has included a 
study of the FSS calibration panel responses. Each of the calibration 
panel runs for each date were. normalized with respect to sun angle 
(see equation 2.10-1) and plotted 

CV^ = V^/cos 0^ 2.10-1 

- FSS response to calibration panel at time t, 

CV^ =5 Sun angle adjusted, FSS response to calibration panel 
at time t, 

6^ « solar zenith angle at time t, 


where 
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See Figures 2,10“‘l and 2,10-2 for examples, Ideally the sun angle 

I 

adjusted response curves would fall on the same curve. This assumes 
that the calibration panel is a Lambertian reflector, which is a good 
assumption down to 50° off normal. The situations which will cause 
the curves to be quite different are; 

1) Data being collected on a cloudy day. 

2) Not filling instrument field of view with the calibration 
panels . 

3) Instrument instabilities . 

4) Shadowing of the calibration panel. 

It should be understood that a review of these curves does not 

necessarily give the whole story about all of the data collected on 

that day. However, a review of the plots does point out particular 

flightlines which the analyst may want to pay closer attention to. 

The values given with the plots (coefficients of variation) attempt 

to place some quantitative value to the variation of the curves so 

that dates can be compared. A summary of the average coefficients of 

variation is given in Table 2,10-1. The average coefficient of variation 

is not computed with band centers 1.425, 1.475, 1.825, 1,875, and 1.925 pm 

since they represent regions of strong water absorption. 

> 

The results in Table 2.10-1 indicate that there is significantly 
less variation on the adjusted calibration responses for the 1975-76 crop 
year than for the 1974-75 crop year. Some of the lower variation can 
be attributed to the shorter span in time between the first and last 
calibration panel observation for the 1975-76 crop year, since data were 
not collected over the extension flightlines. 
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Truck-Mounted Spectrometer Data 

Information, about the repeatability of the spectrometer measurements 
can be obtained using the gray panel data collected at the calibration 
location of the intensive test sites. The reflectance of the gray panels 
does not change appreciably through a year. 

The variation in the measurements represents instrument instabilities, 
data collection and calibration variations, and panel reflectance 
variations. The panel reflectance variations can be due to either 
reflectance changes or differences due to an irregular ground surface (for 
those panels laid on the ground) . 

This information is being compiled and will be included in the 
verification report. Information is lacking, however, for the NASA/JSC 
interferometer (the FSAS system) , The only data available to date is for 
gray panel 1, It is strongly recommended that the FSAS system collect 
measurements over the other four gray panels more frequently than has 
been the case in the past. 

Spectrometer Data Correlation ' 

The spectrometer data is being correlated using the gray panel 

i 

data. A significant part of this data set is planned to be collected 
during mid-May at Finney County, Kansas - the sensor comparison ex- 
periments. The experiments planned will make it possible to compare 
the Purdue/LARS Exotech 20C and Exotech 100 systems and the NASA/JSC 
FSAS and FSS systems at reflectance levels from 5 to 55% with data col- 
lected at the same place and nearly the same time. 
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Table 2.10-1 Summary of FSS Calibration Panel Response Variations 




Williams Co,, 



Finney Co . , Kansas 

. N. Dakota 

Hand Co. 

9 So * 03-lc • 


* 

Ave, Coef^ 


Ave. Coef, 


Ave. Coef. 

Date 

Variation 

Date 

Variation 

Date 

Variation 



1974-75 

Crop Year 



11/05/74 

.13 

6/5/75 

.09 



3/20/75 

.09 

6/22 

.08 



4/8 

,09 

7/10 

.09 



5/14 

.03 

7/18, 

.06 

Not 

5/21 

.12 

7/27. 

.02 



6/2 

.04 

8/5 ^ 

.04 

Applicable 

6/9 

,10 

8/15 

.04 



6/17 

.09 

8/23 

.06 



6/26 

.04 

9/1 

.08 



7/6 

.09 

■ 






1975-76 

Crop Year 



9/16/75 

,05 

5/13/76 

.01 

10/15/75 

.06 

10/3 

.24 

5/28 T 

.01 

10/30 

- 

10/21 

,04 

6/17 

,01 

11/5 

.04 

11/11 

.03 

6/25 

.01 

5/11/76 

.03 

3/18/76 

.02 

7/6 

.03 

6/1 

.03 

3/31 

.05 

7/20 

.03 

6/19 

.05 

4/18 

.02 

7/28 

.03 

7/8 

,04 

5/6 

.03 

8/9 

,02 

7/31 

.03 

6/12 

.06 

8/19 

.08 



6/30 

.07 







1976-77 

Crop Year 



9/28/76 

.01 



9/21/76 

.02 

10/15 

.03 





A 

Average coefficient of 

variation 

(standard deviation/mean) 

for all 

FSS bands 

except those 

at 1.425, 

1.475, 1.825, 

1,875, and 

1.925 ym. 






FIGURE 2.10-1 Z.IO 5 

GRRPH OF SUN RNGLE CORRECTED CRLIBRRTION VRLUES FOR 8/23/75 
SITE WILLIRMS COUNTY N.D. LRTITUDE 48.310N LONGITUDE 103.335W 


SYMBOL 

FLT. LINE 

CRL TIME 

IFOV 

ZENITH 

COSINE ZENITH 

1 

1 

18.21.09 

1.25 

37.3 
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2 

2 

18.39.23 

1.25 

36.7 

0.801 

3 
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19.03. 15 

1.25 

36.6 
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4 

4 

19.22-21 

1.25 

37.0 

0.799 

RVERRGE COEFFICIENT 

OF VRRIRTION (COEF. 

VRR. ) IS 

.06 
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OF CENTER 0.0 0.2 0.4 0.6 0.8 
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0.04 
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0.04 
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0.03 
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GRAPH OF SUN ANGLE CORRECTED^CALIBRATION VALUES FOR 5/6/76 
SITE FINNEY CO. 2. KANSAS LATITUDE 38.148N LONGITUDE 100.665W 


SYMBOL 

FLT. LINE 

CAL TIME 

IFOV 

ZENITH 

COSINE ZENITH 

1 

a 

14.57.51 

1.25 

52.8 

0.604 

2 

9 

15.20.54 

1.25 

48.4 

0.664 

3 

10 

15.42.44 

1.25 

44.2 

0.717 

AVERAGE COEFFICIENT 

OF VARIATION (COEF. 

VAR. ) IS 

.03 


COEFFICIENT BRNO 

OF CENTER 0.0 

VARIATION (MICROMETERS) 


RELATIVE SUN ANGLE CORRECTED RESPONSE 


0.2 


0.4 


0.6 


0.8 


1 


0.02 

0.410 

0.03 

0.430 

0.02 

0.450 

0.02 

0.470 

0.02 

0.490 

0.01 

0.510 

0.01 

0.530 

0.02 

0.550 

0.02 

0.570 

0.02 

0.590 

0.01 

0.610 

0.01 

0.630 

0.01 

0.650 

0.01 

0.670 

0.01 

0.665 

0.00 

0.695 

0.01 

0.710 

0.01 

0.730 

0.01 

0.750 

0.00 

0.770 

0.01 

0.790 

0.01 

0.810 

0.01 

0.830 

0.01 

0.850 

0.01 

0.870 

0.01 

0.890 

0.01 

0.910 

0.01 

0.930 

0.02 

0.950 

0.01 

0.970 

0.02 

0.990 

0.03 

1.010 

0.03 

1.030 

0.05 

1.050 

0.05 

1.070 

0.06 

1.090 

0. 10 

1.125 

0.08 

1,175 

0.06 

1.225 

0.06 

1.275 

0.08 

1.325 

0.20 

1.375 

0.17 

1.425 

0.11 

1.475 

0.06 

1.525 

0.05 

1.575 

0.05 

1,625 

0.05 

1.675 

0.05 

1.725 

0.07 

1 775 

0.10 

1.825 

0.12 

1.875 

0.06 

1.925 

0.07 

1.975 

0.05 

2.025 

0.04 

2.075 

0.04 

2.125 

0.03 

2.175 

0.03 

2,225 

0.03 

2.275 

0.03 

2.325 

0.03 

2.375 
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2.11 Soil Classification and Survey 

Spectral Relationships — Outdoor Exotech Experiment . The Exotech 
Model 20-C spectroradiometer will be used to measure the spectral reflectance 
from surface soils at the Purdue University Agronomy Farm during the week 
of May 9-13, 1977. An experiment lias been set up to measure the effects of 
organic surface residue (com stover) and soil moisture content on the 
reflectance of two soils differing greatly in surface soil color, organic 
matter content, and natural drainage. 

Soils to be investigated are Chalmers silty clay loam, a dark colored 
soil developed under tall prairie grass, and Flncastle silt loam, a light 
colored soil developed under forest vegetation. The Chalmers soil is very 
poorly drained and has a dry Munsell color of 10YR4/1 in the surface horizon, 
and an organic matter content of about six percent. The Fincastle soil is 
moderately well drained and has a surface horizon with a dry Munsell color 
of 10YR6/2 and an organic matter content of about three percent. Chalmers 
soils occur in depressions on broad, gently undulating upland till plains on 
loam to light clay loam glacial till. Fincastle soils occur on broad, 
gently undulating divides on the border of the prairie uplands from a shallow 
covering of silt overlying glacial till. According to the USDA/SCS system 
of Soil Taxonomy (1), Chalmers silty clay loam is classified as a fine loamy, 
mixed, mesic. Aerie Ochraqualf. 

Tbo plot sites have been selected, representing the two soils under 
investigation. On each site twelve plots measuring ten feet by ten feet each 
have been marked off, providing three replications of the treatment combinations. 
Plots have been aligned along the northern edge of the Agronomy Farm bulk 
fields to facilitate spectral measurements at an altitude of about 3 m above 
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the soil surface. 

Treatment combinations consist of two levels of moisture content along 
with two soil surface conditions: with and without organic residue. Soil 

moisture differences will be created by saturating the soil with water one 
day before reflectance measurements are taken on half the plots. The other 
half will remain at the field moisture level. Organic residue of corn stover 
will be applied at a rate of about one-half ton per acre to half the plots. 
This amount of corn residue is about twice that amount necessary to reduce 
erosion (2), while not obscuring the soil background to a large extent. 

The four resulting treatment combinations for each soil are: 

1) moist, bare soil 

2) drier, bare soil 

3) moist soil with corn stover residue 

4) drier soil with corn stover residue 

Surface roughness will be minimized to reduce shadow effects that may 
complicate reflectance measurements. While a“ rough surface promotes better 
infiltration, research indicates that surface residue is the more important 
factor in soil tillage considerations for water conservation (3) . 

Spectral measurements will be made with the Exotech Model 20-C over the 
spectral range from 0.5 to 2.38 ym. Using the 15° FOV mode, the soil surface 

area viewed at a height of 3 m will be circular with an 0.8 m diameter. 

Readings will be taken on clear, cloud-free days as close to solar noon 
as possible to reduce sun angle effects. Measurements will be repeated on 
successive days to verify the repeatability of the spectral responses. 

Surface soil samples of the upper one centimeter of soil collected at 
the time of the Exotech readings will be analyzed for moisture content. 
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organic matter content, and other physical/chemical properties. 

Expected results from this experiment are an increased understanding of 
the effects of soil organic matter, surface crop residue, and soil moisture 
content on spectral reflectance of surface soils. Possible "masking" effects 
of any one of these elements on any of the others will be observed and 
quantified. This should go a long way towards explaining the observed 
spectral differences of similar soils using Landsat data. Cooperative 
efforts with soil survey personnel- will gain from a better understanding of 
the factors affecting soil reflectance in the field situation. The ability 
to identify soil conservation tillage practices from a remote position would 
permit soil managers to evaluate arid monitor the extent of these practices. 
Analysis of Landsat Data 

Geometric correction and temporal overlay of three frames of Landsat 
data from Tippecanoe County, Indiana are being completed. The three frames 
being used are: 

Date Landsat ID 

June 9, 1973 132115593 

April 6, 1975 198715405 


June 20, 1976 


251515411 
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2.12 Improved Analysis Techniques for Multiteraporal Data 


INTRODUCTION 


Ob.lectives . The objectives of this task were to develop and test algo- 
rithms and procedures for analysis of multitemporal data sets; to evaluate 
the utility of classifiers of advanced design (such as perfield classifiers 
and layered classifiers) in contributing to effective multitemporal analy- 
sis; and to document and transfer effective techniques to NASA. 

Approach . The most straightforward way to accomplish multitemporal ana- 
lysis is to extend the "conventional" multivariate classification methodology 
to data sets formed by "stacking" or concatenating the multispectral data vec- 
tors from successive acquisitions over a given site. This presumes, of course, 
that the data sets are geometrically registered. This straightforward approach 
is attractive because of its apparent simplicity, and in some instances it pro- 
duces excellent results, although at some computational cost due to the higher 
dimensionality of the data involved (computation is roughly proportional to 
the square of the total number of channels used for classification). 

However, in addition to the computational cost, there are a number of 
compelling reasons to pursue alternative strategies for multitemporal analysis. 
All are directly or Indirectly related to the relatively limited quantity of 
training data typically available in the face of the increased amount of 
training data actually needed for multitemporal analysis. Thus, the analy- 
sis techniques pursued under this task sought to minimize the need for 
training data while maximizing the performance gains derived from the multi - 
temporal analysis. 

The need for Increased amounts of training data arises in many ways. 

To begin with, the joint effects of natural variability in the ground scene 
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together with differential rates of change in the scene with the passage of 
time can result in an enormous profusion of "spectro-temporal" subclasses. 
Adequate training data must be found for each such subclass, and as the di- 
mensionality of the data increases, the number of training data points needed 
per class also increases. Thus there is a two-fold effect. 

The situation can become worse through the loss of reference data 
("ground truth") due to various factors such as loss of the multispectral 
data at one or more of the acquisitions due to cloud cover or data system mal- 
function; and, in situations involving ground visitation, failure of the obser- 
ver to make or accurately record sufficient observations at all specified sites. 

Registration accuracy is also a factor in training sample availability. 
Training fields, which are often fairly small, must be treated as though they 
were even smaller in multitemporal data because of the increased edge effects 
due to registration errors. 

In summary, although increased data dimensionality plus spectrotemporal 

variability in multitemporal data dictate a need for increased amounts of 

training data, there are many factors at work which tend to reduce the amount 

1 

of training data actually available. Our approach to dealing with this di- 
lemma is to use classification methods which permit "temporal decoupling", 
staged, rather than joint, use of multitemporal data, and thereby substantially 
reduce the demand for training data. The results obtained before termination 
of the project at mid-year had begun to demonstrate the effectiveness of this 
strategy. 

Accomplishments . The accomplishments under this task, a six-month in- 
vestigation terminated by the sponsor at the end of November, 1976, are 
summarized as follows. 

Research was completed on techniques for registration of multitemporal 
data and a comprehensive technical report was prepared and delivered to 
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NASA [1]. Briefly, several registration methods were evaluated relative 
to their geometric accuracy capability. Also an investigation of the 
effects of relative geometric distortion on registration accuracy was con- 
ducted. 

Four multitemporal data sets were assembled. Including registered 
LANDSAT data and comprehensive reference data. All were agriculturally 
oriented, one from the earlier CITARS project and three from LACIE, More 
information on these data sets is contained in the quarterly report for 
the period September 1, 1976 to November 30, 1976. That report also con- 
tains the experimental procedure to have been pursued in evaluating the 
proposed multitemporal analysis- strategies. 

Two methods of multistage classification for multitemporal analysis 
were proposed: the Layered Classifier or decision tree classifier developed 

earlier [2,3], and an extension of the maximum likelihood classifier which 
we have dubbed the Cascade Classifier . The latter is described more fully 
later in this report. 

Due to the premature termination of this task, it was not possible to 
develop or evaluate these approaches comprehensively as had been proposed. 
The preliminary results obtained, however, indicate that both approaches 
have much to offer in the raultitemporal domain. These results were ob- 
tained through a student project in the School of Electrical Engineering 
at Purdue University. 

Conclusions . Despite the father apparent need and potential payoff, a 
comprehensive Investigation of multitemporal analysis strategies has yet to 
be completed. The "stacked vector" method of analyzing multitemporal data 
has definite shortcomings. Multistage methods clearly have promise, based 
on actual analysis results obtained to date. 
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DESCRIPTION OF THE RESEARCH 

Multitemporal Image Registration . Proposed research has been completed 
on techniques for registration of multitemporal multispectral scanner imagery 
(see reference [1]). Several preprocessing techniques were evaluated to de- 
termine what methods provided the maximum improvement in image registration 
accuracy. It was determined that magnitude-of-the-gradient enhancement per- 
formed best. Several correlation measures were also evaluated, and although 
the correlation coefficient gave the best performance, the absolute difference 
method gave almost as good perfoannance with much less computational cost. 

Thus this method is recommended. 

Several other aspects of the registration problem were considered in 
this work. A theoretical method for estimating the variance of the registra- 
tion error was developed. This result was used to derive a correlator which 
optimizes performance. An investigation of the effects of relative geometric 
distortion on registration accuracy was conducted using various signal and 
noise assumptions. Finally the statistical properties of temporal changes in 
LANDSAT data were evaluated. These studies support the advancement of the 
state of the art in image registration and are a part of a growing body of 
knowledge related to the image registration problem. 

The Layered Classifier . Documentation of the basic Layered Classifier 
program, which Implements any user-specified layered decision logic, was es- 
sentially completed. A user-oriented control card summary is available and 
flow charts for the software have been drafted. This LARSYS-compatible pro- 
gram is available to any user of the Purdue/LARS remote terminal system. For 
more information, see Section 2,1 of [3], 


The Cascade Classifier. The Cascade Classifier is a maximum likelihood 
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classifier reformulated and extended to permit "temporal decoupling" or time- 
staged analysis. It is a very versatile classifier model, applicable to mul- 
titype data (e.g., from multiple sources) as well as multitemporal. 

Briefly, the method is derived for the bitemporal case by assuming a 
risk function of the form 

m 

X ^ ^(ijj) p(W2^) jXijXa) 

1 > 2 

where co„. stands for class 1 at time 2. L (j) is the average loss in- 

Zl Ai,A2 

curred by classifying the point characterized by Xi = X(ti) and X 2 = X(t 2 ) 
into class j. The Bayes’ strategy is to minimize L; the maximum likelihood 
strategy develops from taking the loss function X(i,j) to be the familiar 0-1 
loss function. 

Pursuing this development leads to a discriminant function of the form 

J 

d = S p(Xi,X2|ui,U2) pCWzjwi) p(Wl) 

<1)1 

where all subscripts refer to time of acquisition. Notice that p((i)z]a3i) can 
be interpreted as a transition probability. 

By making certain approximations (which are exact under proper temporal 
situations) , this discriminant function can be written as 

d = S p(Xi|l0i) p(X2|w2) p(W2|wi) p(wi) 

0)1 

and it is from this expression that the term "Cascade Classifier" arises. 

This expression implies a classifier model which is shown schematically 
on the following page (Figure 1) . 

There are two advantages to this form of the classifier model, both 
related to "temporal decoupling’’. First, the required density functions 
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are and n 2 -dlmensional, rather than being n^+n 2 -diraensional as in 
the "stacked vector" case. Second, this model allows classification to be 
performed at time 1 and the statistical results passed along to subsequent 
analysis at time 2. Thus optimal use can be made of the data as it becomes 
available, rather than waiting until after the last acquisition. This has 
implications for computational efficiency as well. 

This model is easily generalized to more than two stages or times; 
we simply "cascade" similar stages. 

Software implementing this multistage classifier has been developed. 


Experimental Results . Comparative analyses were performed for the 
Fayette County and Grant County test sites. 

The Fayette County, Illinois, test site was classified into corn, 
soybeans, wheat, and other. The Layered Classifier was designed to use 
t 2 data to resolve uncertainties resulting from classification of t^ data. 
The Cascade Classifier was designed using highly simplified assumptions 
about the prior probabilities and transition probabilities. The results 
were as follows: 


t^ unitemporal (June 29, 1973) 
t 2 unitemporal (July 17, 1973) 
Layered Classifier 

i 

Cascade Classifier 


68% correct overall 
72% correct overall 
82% correct overall 
84% correct overall 


The Grant County, Kansas, test site was classified into alfalfa, corn, 
fallow, pasture, wheat. The same design strategies were used as in the 
Fayette County case. The spectral discrimination problem was clearly much 
more difficult (notice that the; unitemporal performance is considerably 
below that of the Fayette County case) . The results were as follows : 
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unitemporal ( May 9, 1974) 
t 2 unitemporal (July 20, 1974) 
Layered Classifier 
Cascade Classifier 


62% correct overall 
55% correct overall 
68% correct overall 
64% correct overall 


In this case, the Cascade Classifier was apparently much more sensitive to 
the simplified assumptions used for the prior probabilities and transition 
probabilities. 

In neither of the above cases was there sufficient training data avail- 
able to permit multitemporal analysis by the "stacked vector" approach. 


Conclusions and Recommendations . Many questions remain to be answered 
with respect to how to best make use of the information contained in the 
temporal dimension when multitemporal data is available. The "stacked vec- 
tor" analysis approach is the most straightforward to apply, but though 
successful in some Instances, it suffers some severe disadvantages when the 
availability of training data is relatively limited. It is also computa- 
tionally very expensive unless subsets of the available multitemporal chan- 
nels are selected. 

Multistage analysis, using 'methods such as the Layered Classifier and 
the Cascade Classifier, ameliorate many of the difficulties encountered 

I 

with the "stacked vector" approach. The dimensionality of each decision 
i’s minimized, and this reduces both the computation and the amount of 
training data required. Experimental results obtained using these clas- 
sifiers demonstrate that further Investigation of this staged approach to 
multitemporal analysis is warranted. Pursuit of this research is strongly 


recommended. 
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